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Abstract 

Upper limits on the cross-section of the pair-production process e+e^— S'h'^A'^, 
assuming 100 % decays into hadrons, are derived from a new search for the h°A°^ 
hadrons topology, independent of the hadronic flavour of the decay products. Sear- 
ches for the neutral Higgs bosons h° and A°, are used to obtain constraints on 
the Type II Two Higgs Doublet Model (2HDM(II)) with no CP violation in the 
Higgs sector and no additional non Standard Model particles besides the five Higgs 
bosons. The analysis combines LEPl and LEP2 data collected with the OPAL 
detector up to the highest available centre-of-mass energies. The searches are sen- 
sitive to the h°, A'^^qq, gg, r+r^ and h°— ^A^A" decay modes of the Higgs bosons. 
The 2HDM(II) parameter space is explored in a detailed scan. Large regions of 
the 2HDM(II) parameter space are excluded at the 95% CL in the (mh, ^a), (mh, 
tan/3) and (ttia, tan/?) planes, using both direct neutral Higgs boson searches and 
indirect limits derived from Standard Model high precision measurements. The re- 
gion 1 < mh< 55 GeV and 3 < mp^< 63 GeV is excluded at 95 % CL independent 
of the choice of the 2HDM(II) parameters. 
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1 Introduction 



The data collected by the OPAL detector at LEP during the years 1999 and 2000 at 
centre-of- mass energies \/s fsi 192, 196, 200 — 209 GeV are combined with the data 
at the pole, s/s ^ 183 GeV and 189 GeV, to search for neutral Higgs bosons [1-3] 
in the framework of the Type II Two Higgs Doublet Model (2HDM(II)) [4,5] with 
no CP violation in the Higgs sector and no additional particles besides those arising 
from the Higgs mechanism. This study updates the results of a previous OPAL 
pubhcation [6], which included data at the Z° pole, ^/s ^ 183 GeV and 189 GeV. 

In the minimal Standard Model (SM) the Higgs sector comprises only one com- 
plex Higgs doublet [1] resulting in one physical neutral Higgs scalar whose mass is 
a free parameter of the theory. However, it is important to study extended mod- 
els containing more than one physical Higgs boson in the spectrum. In particular. 
Two Higgs Doublet Models (2HDMs) are attractive extensions of the SM since they 
add new phenomena with the fewest new parameters; they satisfy the constraints 
of p ^ 1 [7] and the absence of tree-level flavour changing neutral currents, if the 
Higgs-fermion couplings are appropriately chosen. 

In the context of 2HDMs the Higgs sector comprises five physical Higgs bosons: 
two neutral CP-even scalars, h° and H° (with rrih < mn), one CP-odd scalar, A°, 
and two charged scalars, H='=. The four Higgs masses are free parameters of the 
model. 

Within 2HDMs the choice of the couplings between the Higgs bosons and the 
fermions determines the type of the model considered. In the Type II model the 
first Higgs doublet couples only to down-type fermions and the second Higgs dou- 
blet couples only to up-type fermions. In the Type I model the quarks and leptons 
only couple to the second Higgs doublet. The Higgs sector in the minimal super- 
symmetric extension of the SM [7,8] is a Type II 2HDM, in which the introduction 
of supersymmetry adds new particles and constrains the parameter space of the 
Higgs sector of the model. The 2HDM(II) Higgs potential and detailed description 
of the physical parameters of the model are given in [4,5,7]. 

At the ccntrc-of-mass energies accessed by LEP, the h° and A'' bosons are ex- 
pected to be produced predominantly via two processes: the Higgs-stmhlung pro- 
cess e+e~^h°Z'^ and the pair-production process c+c^^h'^A'^. The cross-sections 
for these two processes, cxhz and cxhA, are related at trcc-lcvcl to the SM Higgs- 
strahlung production cross-section by the following relations [7] : 

e+e-^h^Z" : a^z = sin^ (/?-«) (1) 
e+e-^hOAO : a^p, = cos\(3 - a) X 4|^ (2) 

where af^ is the Higgs- strahlung cross-section for the SM process e+e'^Hgj^jZ^, 
a is the Higgs mixing angle, tan/5 is defined in terms of the ratio of the vacuum 
expectation values, vi and V2, of the two scalar fields, tan/3 = V2/V1 and A = 
A^^^/{A2(f [12m|/s -I- Azh]} accounts for the suppression of the P-wave cross-section 
near the threshold, with Xij — {1—mf /s+nij / s)^ — 4m?m^ / being the two-particle 
phase-space factor. 

In a 2HDM the production cross-sections and Higgs boson decay branching 
ratios are predicted for a given set of model parameters. The coefficients sin^(/3 — a) 
and cos^(/3 — a) which appear in Eqs. (1) and (2) determine the production cross- 
sections. The decay branching ratios to the various final states are also determined 
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by a and In the 2HDM(II) the tree-level couphngs of the h° and A'^ bosons to 
the up- and down-type quarks relative to the couplings of the SM Higgs boson to 
the corresponding fermions are [7] 

,n cos a , n, ^ sin a .n ^ . n, ^ 

h°cc : -, h°bb : -, A°cc : cot /3, A°bb : tan/3, 3 

sm p cos p 

indicating the need for a scan over the range of both angles when considering the 
different production cross-section mechanisms and final state topologies. 

In the analysis described in this paper, detailed scans over broad ranges of 
these parameters are performed. Each of the scanned points is considered as an 
independent scenario within the 2IIDM(II), and results are provided for each point 
in the (mh, mA, tan/3, a) space. The masses mh and ttia are varied such that the 
kinematically accessible range at LEP is fully covered. The choice < /3 < 7r/2 
is derived from V\^V2 > which in the MSSM implies that — 7r/2 < a < [9]. 
This range was studied in [6] to cover an MSSM oriented 2HDM(II). However, 
to completely cover any 2IIDM(II), a has to be varied over an arbitrarily chosen 
angular range of tt. In order to extend the analysis done in [6] beyond the MSSM- 
like 2HDM(II), the domain — 7r/2 < a < n/2 is explored in the present study. The 
final-state topologies of the processes (1) and (2) are determined by the decays of the 
Z°, h° and A° bosons. Higgs bosons couple to fermions with a strength proportional 
to the fermion mass, favouring the decays into pairs of b-quarks and tau leptons 
at LEP energies. However, with values of a and tan/3 close to zero the decays into 
up-type light quarks and gluons through quark loops become dominant, motivating 
the inclusion of flavour independent analyses. 

Section 2 contains a short description of the OPAL detector and the Monte Carlo 
simulations used. A new analysis at the highest LEP energies with improved sensi- 
tivity for the process e+e~— >h°A°— > hadrons, independent of the hadronic flavour 
of the decay products, is presented in Section 3. In Section 4 the confidence level 
calculation method is described. A limit on the cross-section of the pair-production 
process e+e~^h°A*', assuming 100 % decays into hadrons, is given in Section 5. 
The data samples, the final state topologies studied and the external constraints 
used for the 2HDM(II) interpretation are described in Section 6. The 2HDM(II) 
interpretation of the searches is presented in Section 7, and in Section 8 the results 
are summarised and conclusions are drawn. 



2 OPAL detector and Monte Carlo samples 

The OPAL detector [10] has nearly complete solid angle coverage and excellent 
hermeticity. The innermost detector of the central tracking is a high-resolution sih- 
con microstrip vertex detector [11] which lies immediately outside of the beam pipe. 
The silicon microvertex detector is surrounded by a high precision vertex drift cham- 
ber, a large volume jet chamber, and 2;-chambers to measure the z coordinates^ of 
tracks, all in a uniform 0.435 T axial magnetic field. The lead-glass electromagnetic 
calorimeter and the presampler are located outside the magnet coil. It provides, 

^OPAL uses a right-handed coordinate system where the +z direction is along the electron 
beam and where +x points to the centre of the LEP ring. The polar angle, 9, is defined with 
respect to the +z direction and the azimuthal angle, (j), with respect to the horizontal, +x direction. 



Year 


[GeV] 




1999 


191.6 


28.9 


1999 


195.5 


74.8 


1999 


199.5 


77.2 


1999 


201.6 


36.1 


2000 


206.0 


208 



Table 1: The integrated luminosities, £, at each centre-of-mass energy, y^, used 
for the flavour independent h'^A*^— > hadrons search. 

in combination with a number of forward detectors and the silicon-tungsten lumi- 
nometer [12], geometrical acceptance down to 25 mrad from the beam direction. 
The silicon-tungsten luminometer serves to measure the integrated luminosity us- 
ing small angle Bhabha scattering events [13]. The magnet return yoke is instru- 
mented with streamer tubes and thin gap chambers for hadron calorimetry and is 
surrounded by several layers of muon chambers. 

Events are reconstructed from charged particle tracks and energy deposits ( "clus- 
ters" ) in the electromagnetic and hadron calorimeters. The tracks and clusters must 
pass a set of quality requirements similar to those used in previous OPAL Higgs 
boson searches [14]. Charged particle tracks and energy clusters, satisfying these 
requirements, are associated to form "energy flow objects" . A matching algorithm 
is employed to reduce double counting of energy in cases where charged tracks point 
toward electromagnetic clusters [14]. The energy flow objects are then grouped into 
jets and contribute to the total energy and momentum of the event. The association 
into jets is performed by the Durham jet finder algorithm [15]. 

For the h^A"— ^hadrons analysis, the data are separated into four ^/s bins for the 
1999 data, with average ^/s values of approximately 192 GeV, 196 GeV, 200 GeV, 
and 202 GeV. All of the 2000 data is treated together, with an average y/s of 
206 GeV. The luminosities of each of the five data samples are given in Table 1. 

A variety of Monte Carlo samples has been generated in order to estimate the 
detection efficiencies for Higgs boson production and background from SM pro- 
cesses. Monte Carlo signal samples are generated using HZHA [16] on a grid in the 
(mh, ttia) plane, as shown in Figure 1. One thousand signal Monte Carlo events 
are generated and simulated for each grid point for each of the five centre-of-mass 
energies given in Table 1. Samples are generated for four flavour combinations of 
the decays, h^A^^bbbb, h°A°— >-bbcc, h°A°— i^cccc, and h°A°— >• gggg. 

For the background processes the following event generators are used: KK2f [17] 
for (Z/7)*^qq(7), ^^^^""(t) and t+t-(7), BHWIDE [18] for e+e-(7) and KORALW [19] 
and grc4f [20] for four-fermion processes. The KK2f prediction for qq(7) were 
compared to PYTHIA [21] and HERWIG [22] samples. The KORALW prediction for 
hadronic and semi-leptonic four-fermion processes (with no electron in the final 
state) were compared to grc4f samples. JETSET [21] is used as the principal model 
for fragmentation. 

The detector response to the generated particles is simulated in full detail [23] . 



3 Flavour independent search for e+e — )>h^A^ 

For some values of the parameters specifying a 2HDM(II), e+e^^h°Z° is sup- 
pressed, either kinematically or due to small sin^(/3 — a). The decays h*^— i>bb and 
A°— >bb may be also suppressed in a subset of these models because of reduced 
couplings. In such models, the largest signal for Higgs boson production may be 
e"'"e~^h°A*^ where both the h° and the A° decay hadronically, but not necessarily 
to bb. The dominant decay modes of the Higgs bosons may even be to pairs of 
gluons. The final state investigated here is four well-separated jets of hadrons of 
any flavour. 

The search presented here is based on the search procedure applied to the 
189 GeV data [6]. The search is extended by including all data collected at ^/s = 192 
to 209 GcV, and by introducing a likelihood discriminant to combine information 
carried by several different kinematic variables which are measured in each event. 
The results of the search for e+e~— >h°A°— > hadrons in [6] are combined with the 
results of the newer searches; the older data are not re-analysed. 

Without flavour tagging, and without a fixed mass constraint such as the Z'' 
mass, the assignment of dijets in selected candidate events to the h° and A° is 
ambiguous. There are six possible assignments of jets to bosons. The pairing 
chosen in this analysis is the same as that used in [6]. This is the pairing which 
minimizes the of a beam energy and momentum-constrained kinematic fit to the 
(mh, mjC) hypothesis under study. The pairing of each event therefore depends on 
the test masses. 

The Standard Model backgrounds for e"'"e~— ^h^A''— > hadrons are large. One 
of the main sources is e"^e"^qq (approximately 18 pb for events with less than 
20 % of the centre-of-mass energy in initial state radiation at ^/s — 205 GeV [24]), 
and includes events which may have one or more initial-state-radiation photons. 
Much of this background has only two jets. Hard gluon radiation in a fraction of 
e+e"— >qq events produces four-jet final states. In general, the radiation of a gluon 
produces jets close in angle to other jets, and so this background tends to mimic 
e'^e"— >h°A°— > hadrons where either the h° or the A° is light. Conversely if both 
the h" and the A° are assumed to be light, then signal events are more two-jet- like, 
and closely resemble a larger fraction of the qq background. The other main source 
of background is the production of pairs of vector bosons which decay hadronically. 
The most important process is e'^'e"— >W"'"W~— > hadrons (approximately 8 pb at 
^/s — 205 GeV [24]); the process e'^e"— >Z°Z°^*''— > hadrons contributes at a smaller 
level (approximately 0.5 pb at y/s = 205 GeV [24]). These background processes 
produce four-jet events with large invariant masses when jets are combined together 
in pairs. The W^W^ background mimics a signal with mh ~ rriA ~ ttl^. Further- 
more, because the pairing of jets to bosons is ambiguous, the W+W~ and Z°Z°^*^ 
backgrounds can contribute everywhere in the (mh, mA) plane. 

The event selection starts with a cut-based preselection and proceeds with a 
selection based on a hkelihood function. 

3.1 Preselection 

The preselection is based on the preselection used in the search for the Standard 
Model Higgs boson in the Higgs-strahlung process in the four-jet final state [14], 



without the requirement that two jets are consistent with mz- The cuts are: 

1. Events must satisfy the hadronic final-state requirement of [25]. The effective 
centre-of-mass energy, obtained by kinematic fits assuming that initial 
state radiation photons arc lost in the beampipe or seen in the detector [25] 
must be at least 80 % of y^. The value of the jet resolution parameter, 
^34, at which an event is reclassified from 3~ to 4-jet event by the Durham 
algorithm [15] must exceed 0.003. 

2. The C parameter, which gives a measure of the spherical shape of the event [26] , 
must be larger than 0.25. 

3. The probability of a 4-constraint (4C) kinematic fit requiring energy and 
momentum conservation must be greater than 10~^. 

4. The event is forced to have four jets and each of the four jets must have at 
least one charged particle track. 

5. No jet-pairing combination may have a 6-constraint (6C) kinematic fit proba- 
bility greater than 0.2, where the fit constrains the total energy, momentum, 
and the masses of both dijets to mw± ■ 

Table 2 shows the numbers of events passing each of the preselection requirements 
in the data taken in 1999 and 2000, along with the expected backgrounds and the 
efficiency for a signal with — 50 GeV and niA — 100 GeV. 

3.2 Likelihood Selection 

The following nine kinematic variables were considered in the likelihood selection: 

• log(Ax^(mh, ^a)). a full description of the procedure for computing this 
variable is given in [6], where it was used as the only discriminant variable. For 
each event, a 4C kinematic fit is performed, constraining the event energy and 
momentum to the centre-of-mass energy and momentum. The Ax^(mh, ttia) 
value is the additional incurred when constraining one pair of jets to have 
invariant mass and the other to have ttia- The pairing is chosen to minimise 
Ax^(mh, niA)- This is the only variable for which the value depends on the 
test mass combination. 

• I cos 6'thrust I ; obtained from the polar angle, ^thrust; of the thrust axis. 

• The event aplanarity A. 

• log(y34). 

• The jet-angle sum Jg [27]. This variable is the sum of the four smallest dijet 
angles. 

• (-E-max ~ -^min) / V^i the difference between the energy of the highest-energy jet 
and that of the lowest-energy jet, divided by the centre-of-mass energy. 



• The jet-charge-signed cos^w This variable is computed using the jet pairing 
which maximizes the 6C fit probabihty to the W+W~ mass hypothesis. The 
quantity 

^tracks 

i=l 

is computed for each jet j, where ^tracks is the number of tracks in the jet, 
Qi is the charge of the ith track in the jet, and = pi ■ hj where pi is the 
three-momentum of the ith track in the jet and nj is the unit vector pointing 
along the jet axis. If jets j and k are paired together to form a W boson 
candidate and jets I and m are paired to form the other W boson candidate, 
then 



cos 9 



w = — z — z — sig" (Qj + Qk-Qi- Qm) , (5) 

Pi + Pk 



where Pj is the three-momentum of the jet j after the 4C-fit and z is the unit 
vector pointing along the electron beam axis. 

• log(VFcco3); the logarithm of the WW matrix element calculated by the EXCA- 
LIBUR program [28] using the CC03 set of diagrams. The four- vectors of the 
jets after the 4C-fit are used as inputs to the calculation. The matrix element 
is computed for all possible assignments of jet pairs to W bosons and the 
largest value is used. 

• log(l^QCD), the logarithm of the e+e~— >qq— > four-jet matrix element [29]. 
The matrix element is computed for all possible permutations of jets and the 
largest value is used. 

Reference histograms are formed for each likelihood input variable, for each 
signal grid point, separately for the e"'"e~— >qq background (2f), the e"'"e~^ qqqq 
background (4f), and the expected signal, accumulating events which pass the pres- 
election requirements. The £^£^qq and iuiqq backgrounds are expected to be small 
after the preselection - their numbers are included in the 4f background numbers 
in Table 2. They do not contribute to the reference histograms, but are accounted 
for in the final background estimates in the likelihood output histograms. The 
backgrounds from two-photon processes are negligible after the preselection. 

Figure 2 shows the distributions of the likelihood input variables for the data col- 
lected in 1999 and 2000, the corresponding background estimate, and the expected 
signal for a fully gluonic decay for the hypothesis mh = 50 GeV, ttia — 60 GeV. 

For any point in the (mh, ttia) plane within kinematic reach of the LEP beam 
energy, and with > 30 GeV and mx > 30 GeV, a separate likelihood function 
may be constructed from the reference histograms of the input variables. These 
are formed by interpolating the signal reference histograms using nearby Monte 
Garlo grid points. The background reference histograms must be interpolated also 
for the log(Ax^(mh, mA)) variable. These interpolations make use of the method 
described in [30], extended to interpolate histograms which are functions of two 
variables, mh and mA- The likelihood output histograms are also interpolated, 
separately for the signal and each background contribution, but each bin's contents 
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2f l)kg 


if bkg 


tot l)kg 


(lata 


off [<X] 


(1) 


1031.6 


3206.4 


4238.0 


4479 


76.6 


(2) 


976.0 


3203.1 


4179.1 


4436 


76.6 


(3) 


887.3 


2977.2 


3864.5 


3807 


71.0 


(4) 


784.7 


2820.2 


3605.0 


3580 


70.4 


(5) 


754.3 


2035.7 


2789.9 


2835 


68.5 



Table 2: The numbers of events passing each of the preselection requirements 
in the data taken in 1999 and 2000, along with the expected backgrounds and 
the luminosity-weighted average efficiency for a signal with mh = 50 GeV, 
ruA = 100 GeV. 

is linearly interpolated between the Monte Carlo grid points. The interpolated 
reference histograms allow the computation of the value of the likelihood function 
for each candidate for each point in the (mh, mA) plane, and the interpolated 
signal and background likelihood histograms allow comparison of the data likelihood 
distribution with the signal and background predictions. 

Signal events with h'^A^— >bbbb are easier to separate from the backgroTind than 
signal events with h°A°— gggg since the gggg case has a poorer reconstructed mass 
resolution which deteriorates the log(Ax^(mh, mA)) likelihood variable. The signal 
reference histograms are created using samples of h°A°— >bbbb signal Monte Carlo 
events, and the signal hkelihood histograms are filled with h'^A'^— > gggg events, 
ensuring the statistical independence of the reference and likelihood histograms, 
and also ensuring the conservativeness of the performance over the possible final 
states of the h° and A° decays. 

Example likelihood distributions are shown in Figure 3 for all data collected in 
1999 and 2000, along with the SM background expectations and signal expectations, 
for three test-mass hypotheses, (mh, mA) = (50 GcV, 100 GeV), (50 GeV, 60 GeV), 
and (30 GeV, 60 GeV). The corresponding distributions of log(Ax^(mh, mA)) are 
also shown to illustrate how the distributions of the signal, the expected back- 
grounds, and the candidates change with the test-mass hypothesis. 

The distribution of the likehhood is used directly as the input to the limit 
calculation. In the presence of systematic uncertainties on the background rate, 
including bins with low expected signal-to-background ratios reduces the sensitivity 
of the search. A lower cut on the likelihood variable of 0.8 is chosen, independent of 
the test mass hypothesis, in order to improve the sensitivity of the search. Table 3 
lists the numbers of events passing the hkelihood cut for each of the test masses on 
the Monte Carlo grid, the expected backgrounds, and the expected signal efficiencies 
for the h^A''^ gggg decay hypothesis. The signal efficiencies are also calculated 
separately for the cccc, bbbb, and the bbcc decay hypotheses. For nearly all test 
mass combinations, the h^A"— gggg hypothesis yields the least efficiency, and for 
the remainder, the differences are within the uncertainties. 
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30.0 
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43 
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40.0 
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50.0 
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32.3 
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70.4±6.6 
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34.0±2.9 
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36.5 


33.6 


80.0 


50.0 


11.9 


57.2 


69.1±6.5 
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32.2 
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100.0 


50.0 


12.6 
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79 
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25.5 


26.0 
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50.0 


18.5 
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78.4ib7.4 
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20.4 


22.7 


60.0 


60.0 


15.2 


43.0 
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40.0 
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16.8 
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70.0 
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36.0 


31.5 


34.6 


80.0 


70.0 


8.7 


99.2 


107.9±10.1 


106 


18.1±1.8 


27.2 


29.8 


28.0 


100.0 


70.0 


12.0 


84.1 


96.1±9.0 


100 


16.4±1.7 


23.3 


24.9 


21.2 


80.0 


80.0 


1.7 


56.9 


58.6±5.5 


54 


2.5±0.5 


4.6 


9.6 


6.1 


90.0 


80.0 


4.4 


92.9 


97.2±9.1 


110 


11.5±1.3 


18.2 


19.5 


14.6 


90.0 


90.0 


21.0 


142.0 


162.9±15.3 


162 


25.7±2.3 


31.7 


32.5 


29.8 



Table 3: The expected Standard Model backgrounds, observed data counts, and 
expected signal efficiencies for the flavour independent h°A° search as a function of 
the test mass hypotheses. All data collected in 1999 and 2000 are combined, and 

the signal efficiencies are luminosity-weighted averages. The efficiencies are listed 
separately for gggg, cccc, bbbb and bbcc decay hypotheses. The errors given in 
table are statistical and systematic added in quadrature. 

3.3 Systematic Uncertainties 

Because of the cut on the likelihood of 0.8, the systematic uncertainties are evaluated 
only for the numbers of events passing this likelihood cut. Correlations between the 
uncertainties on the signal and background rates, as well as between samples taken 
at different centre-of-mass energies are evaluated and used in the computation of 
the confidence levels. 

• Monte Ceirlo Statistics For a typical point in the Monte Carlo test mass 
grid, the MC statistical error on the signal rate is 7 %, but it is larger for 
models for which both mh and are small, due to the lower signal efficiency 

1 1 



for such models. For nearly all model points, the MC statistical uncertainty 
on the background rate is between 1.5 % and 3 %, but for signals when either 
mh or mA is low, it can be as large as 5 %. An overall uncertainty of 7 % and 
3 % is taken for the signal and background rates, respectively. 

• Jet Energy Resolution The jet energy resolution is uncertain in the barrel 
region by about 3-5 %, but this uncertainty is approximately 15 % in the 
endcap regions. Uncertain jet energy resolution results in approximately a 
5 % uncertainty on both the signal and background rates. 

• Jet Energy Scale The jet energy scale is uncertain at the 1 % level. The 
corresponding uncertainty on the rate of events passing the likelihood cuts is 
approximately 2 % for both the signal and the background. 

• Jet Angle Resolution A 22 mrad angle resolution uncertainty in both 6 
and results in an uncertainty on the rate of events passing the likelihood 
cuts that is approximately 1 % for both the signal and the background. 

• Interpolation Uncertainty To test the reliability of the interpolation pro- 
cedure, a single Monte Carlo signal point on the grid is deleted, and the 
interpolation procedure is used to replace it, and the process is repeated for 
all MC grid points not on the edges. Using the differences found in the se- 
lection rates, an uncertainty of 3.5 % is assigned to the signal efficiency and 
an uncertainty of 4 % is assigned to the background rate due to interpolation 
errors. 

• Four-Fermion Cross-Section Nearly all of the background at model points 
near the expected limit is from four-fermion production. It is dominated by 
e+e~— ^W+W" production, although e+e"— >Z°Z°^*^ contributes as well. A 2 % 
uncertainty is assessed on the production cross-section of these events [31]. 

• Monte Carlo Background Sample Comparison The KK2F Monte Carlo 
generator using PYTHIA as the fragmentation and hadronisation model was 
used to generate the central values for the qq background rates. These rates 
were compared to the same KK2F sample but rc-hadronised with HERWIG, and 
also with a sample generated entirely with PYTHIA. KORALW was used to gen- 
erate the qqqq, the qq£~^£^ and the qqiui background rate central values, and 
grc4f was used as the comparison generator. An uncertainty of 5.4 % is 
assessed on the background rates passing the likehhood cuts. 

When all uncertainties are added in quadrature, the systematic error on the 
background rate is 9.4 % and the systematic error on the signal rate is typically 
9.6 %, with larger values for the signal systematic uncertainty for models with low 
efficiency, due to the Monte Carlo statistical uncertainty. The correlations between 
the signal and background uncertainties and between years are taken into account. 
The Monte Carlo statistical error affects the signal and background predictions 
and are uncorrelated between energies, and signal and background. The other 
uncertainties are correlated between centre-of-mass energies. 
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4 Confidence level calculation 



Following the statistical method described in [32,33], the direct searches listed in 
Section 6.1 are combined to increase the Higgs boson discovery potential and, in 
case of absence of signal, the exclusion power. 

The confidence levels are derived from a test statistic, Q, which is defined such 
that Q quantifies the compatibility of the data with two hypotheses: a) the back- 
ground hypothesis, and b) the signal+background hypothesis. The confidence lev- 
els are computed from a comparison of the observed test statistic and its prob- 
ability distributions for a large number of simulated experiments for these two 
hypotheses. The results of the different search channels are expressed in bins of 
discriminating variables defined in the individual searches {e.g. mass, likelihood, 
neural network output, etc.). The ratio Q = Cg+h/ of the binned likelihoods 
for the two hypotheses is chosen as the test statistic. The confidence level for 
the background hypothesis, CLb, is defined as the probability to obtain values of 
Q no larger than the observed value Qobs, given a large number of hypothetical 
experiments with background processes only, CLb = P(Q < Qobs| background). 
Similarly, the confidence level for the signal-|-background hypothesis, CLg+b, is de- 
fined as the probability to obtain values of Q not larger than observed, given a 
large number of hypothetical experiments with signal and background processes, 
CLg+b = P{Q < Qobs I signal + background). In principle, CLg+b could be used 
to exclude the signal+background hypothesis, given a model for Higgs boson pro- 
duction. However, this procedure may lead to cases when a downward fiuctuation 
of the background would allow hypotheses to be excluded for which the experi- 
ment has no sensitivity due to the small expected signal rate. Therefore the ratio 
CLs = CLs+b/CLb is used. It is always greater than CLg+b and the limit obtained 
in this way is thus conservative. We adopt this quantity for setting exclusion lim- 
its and consider a hypothesis to be excluded at the 95 % confidence level if the 
corresponding value of CLg is less than 0.05. 

The expected confidence levels are obtained by replacing the observed data with 
a large number of simulated events for the background only or signal-|-background 
hypotheses. 

The effect of systematic uncertainties for the individual channels is calculated 
using a Monte Carlo technique. The signal and background estimations are varied 
within the bounds of the systematic uncertainties, assuming Gaussian distributions 
of the uncertainties. Correlations are taken into account. These variations arc con- 
voluted with the Poisson statistical variations of the assumed signal and background 
rates in the confidence level calculation. The effect of systematic uncertainties on 
the exclusion limits generally turns out to be small. 

5 Model independent interpretation 

Since no excess of data has been observed, the flavour independent search for 
e"'"e~— >^h°A°, described in Section 3, is used to set 95 % CL upper limits on the 
Y^opo production cross-section assuming 100 % hadronic branching ratios. The 
cross-section for e"'"e~— >h°A° is determined by mh, rriA, and the scale factor c^, 
analogous to the cos^(/3 — a) factor of 2HDMs. The scale factor is defined as the 
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ratio of the pair-production cross-section in the model considered and the pair- 
production cross-section given in Eq. (2) with cos^(/9 — a) = 1. The couphng hmit 
is calculated by finding the value of for which CLg = 0.05, assuming both the h° 
and the A° to decay 100 % hadronically. Figures 5 and 6 show the 95% median 
expected and observed upper limit on c^, respectively, as a function of the test- mass 
hypotheses, {m^, m^). 

The search for e"'"e~^h*'A*'^hadrons is sensitive in the region where the produc- 
tion cross-section exceeds 200 ^ for between 189 GeV and 206 GeV. For 100% 
decays to hadrons and = 1, this corresponds to mh -l- ~ 130 GeV; for these 
test masses, the production cross-section does not change significantly between the 
centre-of-mass energies of the data used. Monte Carlo signal samples were gener- 
ated with mh-l-mA up to 170 GeV, but due to the low signal cross-sections for high 
mass Higgs bosons, the 1 — CLb results are reported only up to mi^ + mx = 145 GeV. 
The higher-mass Monte Carlo samples are needed to interpolate the reference and 
likelihood histograms over the entire range over which the results are produced. 

The composition of the background which passes the hkelihood cut depends 
strongly on the test-masses mh and uia- For low mh and mA, the 2f background 
dominates, but for models near the limit, the 4f background is the most important. 
The signal efficiency after the likelihood cut also depends strongly on the mass 
hypotheses. For mh ~ mA ~ 30 GeV, the 2f background is quite large and closely 
mimics the signal. Only a very small fraction of the signal and of the background 
passes the likelihood selection requirement because of the reduced separation power 
between the signal and the background - the signal efficiency for this particular 
model is only 1.9 ±0.5%. This small efficiency is compensated by the large expected 
signal cross section for low Higgs masses. Another mass hypothesis where the signal 
efficiency is low is mh ~ mA ~ 80 GeV. In this case, the W+W~ background is 
dominant, and the separation between the signal and the background is poor. For 
this model, the signal cross-section is between 70 fb at = 189 GeV and 100 fb at 
y/s = 206 GeV and is therefore beyond the range of sensitivity of the search. Over 
much of the range of tested mass hypotheses, the signal efficiency is between 20% 
and 30%, and signal cross-sections as low as 200 fb are excluded. 

To determine whether a signal was observed, wc compute 1 — CLb for each 
point in the (mh, mA) plane in the search region. Given the mass resolution of 
approximately 3 GeV for the sum of mh and mA, and approximately 7 GeV for the 
difference between mh and mA, there are approximately 160 independent searches 
each of which may have an excess or deficit, diluting the significance of 1 — CLb. 
Thus, more than one independent excess in 1 — CLb is expected at the percent level. 
Figure 4 shows 1 — CLb on a logarithmic scale as a function of mh and mA. Nowhere 
is 1 - CLb below 1%. 

As a result of the improved sensitivity of the analysis and the inclusion of data 
taken at higher centre-of-mass energies in the years 1999 and 2000, the excluded 
domains in Figure 6 are extended substantially beyond those obtained in [6] . 



6 Search channels and external constraints used 
in the 2HDM(II) interpretation 

6.1 Data samples and final state topologies studied 

The present study relies on the data collected by OPAL at ^/s ^ mz and from 
^/s 183 GeV to 209 GeV, the highest e+e" collision energy attained at LEP. This 
paper uses existing published analyses for all but the new h°A*^ flavour independent 
channel described in Section 3. Channels that use b-tagging provide useful infor- 
mation in regions of the 2HDM(II) parameter space where the Higgs bosons are 
expected to decay predominantly into bb pairs. Flavour independent channels, not 
using any b-tagging information, arc also included in the combination in order to 
explore the regions at low a or low tan/5, where the decays of the h*^ and A" bosons 
into bb and r"'"r~ pairs are suppressed. In Table 4 the references to the published 
OPAL papers for the direct search channels combined in the present 2HDM(II) in- 
terpretation are given, together with the corresponding centre-of-mass energies at 
which they were performed. 

The channels used at ^/s ~ mz, 183 and 189 GeV are the same as in [6]. 
The integrated luminosities, the numbers of candidate events, the expected SM 
backgrounds and the efficiencies for each of the b-tagging (flavour independent) 
h^Z^ channels at 192 < < 209 GeV are given in Table 5 (Table 6). 

The detection efficiencies quoted in Tables 5 and 6 are given as examples for 
specific values of mh. 

The integrated luminosities, the numbers of candidate events, the expected SM 
backgrounds and the efficiencies for the most relevant b-tagged h°A° channels are 
given in Table 7. 

When scanning the parameter space the efficiency is calculated for each point 
in the (mh, niA) plane for each of the final states considered. 

6.2 External constraints 

In addition to the combination of the direct search channels, the following external 
constraints are applied in every parameter space point considered: 

(a) A powerful experimental constraint on extensions of the SM is the determina- 
tion of the total width of the Z° boson, Pz, at LEP [37]. Any possible excess 
width obtained when subtracting the predicted SM width from the measured 
value of Pz, APz, can be used to place upper limits on the cross-section of 
Z° decaying, as in the 2HDM, into final states with h° and A° bosons [38]. 
The maximum additional contribution to the total Z° width that is compatible 
with the measured width at 95 % CL is APz = 6.5 MeV, obtained from the 
latest LEP combined Z° hneshape results [24]. An expected increase of the 
partial width of the Z° is evaluated for each scanned parameter space point 
in the 2HDM(II); if it is found to exceed the experimental limit, the point is 
excluded. 

(b) The decay mode independent search for e"'"e~^SZ° [39], where S is any scalar 
particle produced in association with the Z" boson, provides an upper limit on 
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Decay topologies 


mz 


183 


189 


192 ^ 209 




bbqq, hhui/, bbe''"e~, bb/u'^/x", bbr'^'r" 


NA 


[6] 


[6] 


[14] 


qqvi', qqe'^'e", qqfj,'^fj,~, qqr"^r~ 


[6] 


NA 


[6] 


[34] 


r+r~ qq 


[6] 


[6] 


[6] 


[34] 


qqqq 


XA 


XA 


[fi] 


[34] 




qqr'^r or r'^r qq 


[6] 


NA 


NA 


NA 


bbbb, bbr'^r~or r"*'r~bb 


NA 


[6] 


[6] 


[35] 


qqqq 


NA 


NA 


[6] 


This paper 


mh > 2mA, hOAO^A^A^A^ 


bbbl)l)l) 


[6] 


[6] 




[35] 


mh> 2mA and mA < 2mb, h'^ZO^AOAOZ'' 


A°— >-cc, r"^r , gg and Z^^vu, ji^ ^ , e+e 


NA 


NA 


[36] 


[36] 


mh > 2mA, hOZ°^AOA°Z*^ 


bbbbqq, bbbbi/z/ 


NA 


[6] 


[6] 


[35] 


bbbbe+e~, bbbb/i^/i", bbbbr+r" 


NA 


[6] 


[6] 


NA 



Table 4: Direct search channels combined in the present interpretation of the 
2HDM(II). The searches with qq final states include gg production as well. The 
numbers in the table give the references to the OPAL publications where a full 
description of the channel can be found. Channels marked NA do not exist. The 
h^'A^^qqqq (gggg, ggqq) analysis for the data taken in the years 1999 and 2000 is 
new and is described in Section 3 of this paper. 



the scaling factor A;, defined as crsz = ka^, where crsz is the production cross- 
section for a scalar S in association with a Z°, and crf^ is the expected SM cross- 
section for ms = ^HsM- This translates into a limit on the production cross- 
section in each parameter space point of the 2HDM(II) for which Uhz > 
at 95 % CL. 

(c) In regions of the 2HDM(II) parameter space for which 4 < mh (mA) < 12 GcV 
a special study was performed in [40], and 95 % CL limits were obtained on 
the Yukawa couplings of h*^ and A° to down-type fermions. These limits are 
apphed as an external constraint in the scan of the the 2HDM(II) parameter 
space described in Section 7. 

The production of any neutral low mass scalar particle in association with the Z° 
was investigated in [41] and, for mh < 9.5 GeV, a mass-dependent upper limit on the 
Higgs boson production cross-section was obtained. This limit was translated in [6] 
into an upper limit on the production cross-section for mh below 9.5 GeV, which was 
considered as an external constraint in combination with the Z°-width constraint 
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Channel h^ZO^ 


£ [pb-i] 


Data 


Total bkg 


eff [%] 


192 < V« < 202 GeV 


bbqq {rrih = 100 GeV) 


217.0 


30 


28.8 ±4.2 


42.0 




212.7 


10 


13.9 ± 1.6 


46.9 


bbr+r^/r"'"r~qq 


213.6 


5 


5.1 ±0.8 


25.8 


bbc^e" 


214.1 


3 


4.1 ±0.7 


57.2 


bb/i+/Li^ 


213.6 


6 


3.3 ±0.5 


62.5 


200 < V« < 209 GeV 


bbqq (rrih = 115 GeV) 


207.3 


20 


17.5 ±2.6 


40.0 


bbz/P 


207.2 


11 


8.9 ± 1.0 


40.7 


bbr+r~/r^r^qq 


203.6 


5 


4.5 ±0.7 


25.6 


bbe+e~ 


203.6 


1 


3.6 ±0.7 


52.9 


bb//+/i~ 


203.6 


4 


3.4 ±0.5 


59.2 



Table 5: The h°Z° b-tagging channels for data collected in the year 1999 with 
192 < < 202 GeV, and for data collected in the year 2000 with 200 < < 209 
GeV, respectively. The integrated luminosities (£), the numbers of events after the 
final likelihood or Neural Network cut for the data and the expected background, 
normalised to the data luminosity are shown. The errors on the total background 
include modeling uncertainties and Monte Carlo statistics. The last column shows 
the detection efficiencies for a Higgs boson with mh = 100 GeV for the year 1999 
data and = 115 GeV for the year 2000. Since the four-jet channel relies on a 
mass- dependent analysis, the numbers quoted in the table are given as an example 
for mh= 90 GeV (100) for the year 1999 (2000) data. For the four-jet channel, 
the efficiency is computed only for h°— > bb decays, while for the missing-energy, 
electron and muon channels the efficiency is for all decays of the h°, assuming SM 
branching fractions. For the tau channel, the efficiency is quoted for the processes 
Z'^h^— >-T"'"T~(h°— >-all) or Z^h*^— >qqr"''r~ assuming SM branching fractions. 
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Channel h^Z^^ 


JC [pb-^] 


Data 


Total bkg. 


eff [%] 


192 < < 202 GeV 


qqqq (mh = 90 GeV) 


217.0 


290 


290.6 ±37.7 


52.0 


qqui' 


212.7 


68 


70.8 ± 10.8 


45.2 


qqr'^T'' / r'*"r~qq 


213.7 


1 


4.6 ±0.8 


27.1 


qqe+e~" 


214.1 


13 


8.3 ±2.5 


58.6 


qq/i+/x~ 


213.6 


7 


7.6 ±1.5 


64.8 


200 < ^/i < 209 GeV 


qqqq (mh = 100 GeV) 


207.3 


263 


235.3 ±28.2 


55.0 


qqz^P 


208.2 


55 


62.3 ±2.2 


48.4 


qqr^r^/r^r^qq 


205.3 


2 


4.2 ±0.8 


21.0 


qqe+e~ 


208.2 


10 


8.3 ±2.5 


58.8 


qq/i+^- 


207.8 


9 


7.4 ±1.4 


62.4 



Table 6: The h^Z" flavour independent channels for data collected in 1999 with 
192 < < 202 GeV and in year 2000 with 200 < < 209 GeV: the integrated 
luminosities (£), the numbers of events after the final likehhood or Neural Network 
cut for the data and the expected background, normalised to the data luminosity. 
The errors on the total background include modeling uncertainties and Monte Carlo 
statistical errors. The last column shows the detection efficiency for a Higgs boson 
decaying to quark or gluon pairs with mh = 90 GeV and with mh = 100 GeV for 
data collected in year 1999 and 2000, respectively. Since the four-jet channel relies 
on a mass-dependent analysis, the numbers quoted in the table are given as an 
example for mh = 90 GeV (100) for the year 1999 (2000) data. 



Channel h^A^^ 


(mh, mA) [GeV] 


L [pb-i] 


Data 


Total bkg 


eff [%] 


192 < < 209 GeV 


bbbb 


(90, 90) 


424.3 


22 


19.9 ±2.01 


49.4 


bbr'*"r~ / T'*"r~bb 


(90, 90) 


417.2 


13 


13.2 ±2.02 


42.5 


bbbbbb 


(100, 40) 


424.3 


22 


19.9 ± 2.01 


59.4 


199 < < 209 GeV 


bbbbz^z^ 


(100, 40) 


207.2 


19 


17.2 ± 2.30 


66.0 


bbbbqq 


(100, 40) 


207.3 


20 


17.5 ± 2.6 


31.5 



Table 7: The h'^A" channels for data collected in the years 1999 and 2000. The 
integrated luminosities (£), the numbers of events after the final likelihood or Neural 
Network cut for the data and the expected background, normalised to the data 
luminosity. The errors on the total background include modeling uncertainties and 
Monte Carlo statistics. The last column shows the detection efficiency for the mass 
combination given in the first two columns. 
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(item (a) in the previous list). The decay mode independent results included as an 
external constraint (item (b) in the previous list) provide better exclusion power 
than the upper limit on Higgs boson production cross-section obtained in [41], 
increasing the exclusion power in the 2HDM(II) parameter space for low mh values. 
The external constraint on the Yukawa couplings extends the exclusion power to 
regions of the 2HDM(II) parameter space with low values of mh and uia and large 
tan/5. In general, the external constraints applied in the present study improve 
significantly the results that were obtained in [6]. 

7 Two Higgs Doublet Model interpretations 

The interpretation of the searches for the neutral Higgs bosons in the 2HDM(II) 
is done by scanning the parameter space of the model. Every (mh, mA, tan/3, a) 
point determines the production cross-section and the branching ratios to different 
final states. The 2HDM(II) parameter space covered by the present study is: 

• 1 < mh < 130 GeV, in steps of 1 GeV 

• 3 < mA < 220 GeV, in steps of 1 GeV; 
300 < mA < 500 GeV, in steps of 100 GeV; 
0.5 < mA < 2.0 TeV, in steps of 0.5 TeV 

• 0.4 < tan/3 < 40, in steps of 1° in (3, from /3 = 22° to /3 = 88° and an 
additional point corresponding to tan/3 = 40 

• Q;=-7r/2, -7r/4, 0, 7r/4, and 7r/2 

The values of a are chosen to extend the analysis to the particular cases of 
maximal and minimal mixing in the neutral CP-cvcn sector of the 2HDM(II) {a — 
±7r/4 and ±7r/2, respectively) and of BR(h°— s>bb) = (a = 0). The extreme 
cases a = ±7r/2 are equivalent in the 2HDM(II) since the mass matrix of the CP- 
even neutral Higgs sector, containing the Higgs doublets, becomes diagonal. In [6] 
models with 1 < mh < 100 GeV (in steps of 1 GeV), 5 < mA < 100 GeV (in 
steps of 1 GeV, and larger steps up to 2 TeV), 0.4 < tan/3 < 58 (in steps of 1° in 
/3) and a = — 7r/2, — 37r/8, — 7r/4, — n/S and 0, were considered. The present 
study is extended to cover models with positive values of a {a = 7r/4 and 7r/2), 
which are not allowed in the MSSM-like scenarios. Furthermore larger values of mh 
and mA are explored in this study, due to the increased sensitivity of the search 
channels to high mh and mA values: the data analysed have been collected at larger 
centre-of-mass energies and with larger luminosities than those in [6]. For tan/3 < 
0.4, the theoretical predictions become unreliable. For tan/3 > 40 the width of the 
A° and h° becomes non-neghgible. The decay mode independent study introduced 
in Section 6.2. [b] is providing exclusion from mh ~ 1 GeV down to 1 KeV, where 
the limit on the cross-section scaling factor, k, with respect to the Standard Model 
Higgs-strahhmg cross-section is of the order of 0.06 [39]. Below mA ~ 3 GeV 
radiative corrections become unstable inducing large fluctuations in the calculated 
cross-sections. 

In the present study the other two free parameters of the model, mn and mH± , 
are set at values of mn = 210 GeV and mH± = 1 TeV, above the kinematically 
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accessible region at LEP. A scan over values of the masses mn and mH± up to 2 
TeV has been performed and no change has been observed in the production cross- 
sections and branching ratios to final state topologies of h° and A°, as expected 
from the theory. The HZHA Monte Carlo generator [16] that includes the 2HDM(11) 
production cross-sections and branching ratios of the Higgs particles has been used 
to scan the parameter space. This generator includes next-to-next-to-leading order 
QCD corrections and next-to-leading order electroweak corrections. The branching 
ratios and cross-sections obtained were cross-checked with the results of another 
generator [42] in which QCD corrections are computed only up to next-to-leading 
order. The comparison showed an agreement better than 1 % between the results 
of the two programs. 

The results of all the individual search channels at the studied centre-of-mass 
energies are combined statistically to provide 95 % confidence level (CL) limits, 
which are extracted using the method explained in Section 4. By applying the 
external constraints discussed in Section 6.2 additional regions of the parameter 
space are excluded at 95 % CL. Although the flavour independent channels supply 
a unique way of investigating parameter space regions where the branching ratio 
h°-^>bb or A°^bb is highly suppressed {e.g., low a and tan/? regions), they have 
a poor sensitivity with respect to the b-tagging channels outside these regions. 
The use of b-tagging information substantially reduces the background coming 
from W+W~ events and improves the sensitivity of observing Higgs bosons even 
in regions of the 2HDM(II) parameter space where only small branching ratios 
for h°^bb are expected. The fiavour independent and b-tagging searches have 
candidate events in common, as is also true for the b-tag 4-jets h^A^ and h°Z° 
searches. To avoid double-counting of candidate events, in each of the parameter 
space points only the channel that provides the better expected confidence level 
among the ones that have some candidates in common is used for the extraction of 
the limits. 

The direct searches for the process e+e^— >-h°Z° (e+e^^h''A°) in the Z° data 
contribute mainly in the rrih < 50 GeV {rrih < 60 GeV) region. For ^/s > 189 GeV, 
since the flavour independent h°Z° and h°A° analyses have been performed in the 
mass regions rrih > 60 GeV and nih, tua > 30 GeV, respectively, only b-tagging 
channels are applied below these masses. The fiavour independent analyses provide 
exclusion for the whole tan/3 range and for the tan/3 < 1 regions for a = and 
a — ±7r/4, respectively. 

In Figures 7(a-d) the excluded regions in the (mh, ttia) plane are shown for 
the chosen values of a, together with the calculated expected exclusion limits. A 
particular (mh, mA, ct) point is excluded at 95 % CL if it is excluded for all scanned 
values of tan/3. Different domains of tan/3 are studied and described below: a) 
0.4 < tan/3 < 40 and b) 0.4 < tan/3 < 1.0 or 1.0 < tan/3 < 40, for which enlarged 
excluded regions are obtained. 

a) 0.4 < tan/3 < 40 (darker grey area): 

• For all a values, in the region below m-^ < 10 GeV the sensitivity of the chan- 
nels at ^/s ~ mz is poor. The exclusion in this region is due to the application 
of the external constraints, as explained in section 6.2, in particular to the Z 
width constraint. Both the h°Z° and h°A° production processes contribute to 
the natural width of the Z°. The exclusion provided by the h'^Z^ process is 
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valid for any value of mA- On the other hand, the exclusion provided by the 

h°A° process is kinematically limited to the regions where mA + < rriz- 
The contribution of the h°Z° production cross-section to the Z° width depends 
on the argument (/3 — a), and it becomes large enough for this process alone 
to provide exclusion in different tan/3 domains for the a values considered. 

• The shape of the exclusion plot in Figures 7(a) and (b) for mh < 30 GcV is re- 
lated to the kinematic constraint on the h^A^ production in the data, which 
for /3 — q; or TT is the only allowed process, since the h°Z° production cross- 
section vanishes. Since the a values are such that the condition /3 — a ~ or 
TT is never achieved in Figures 7(c) and (d), the domains with mh + mx > \/s 
can be excluded by the h'^Z'^ channels. For mh > 20-30 GeV, the high energy 
data open a new kinematic region and are able to exclude large (mh, mA) ar- 
eas, as can be seen in Figures 7(a) and (b). In the same figures the exclusion 
in the observed rectangular contour 20 < mh < 30 GeV and 90 < mA HO 
GeV is due to the recent optimisation of the e+e^^h°A'^^bbbb analysis in 
the same kinematical region [35]. The unexchided region 20 < mh < 30 
GeV and 60 < mA < 90 GeV is due to a small excess in the data with 
respect to the expected background of about the same size as the expected 
signal, for tan/3 fa 8 in Figure 7(a) and for tan/3 fs 0.7 in Figure 7(b), 
respectively. 

• The (mh, mA) points below the semi-diagonal defined by mh > 2mA, for 
which the process h°— >A°A° is kinematically allowed, can only be excluded 
by the high energy channels for restricted tan/3 ranges depending on the a 
values examined. For very low values of tan/3 the branching ratio for A°— ^bb 
vanishes, causing unexcluded regions for all values of a. However, these are 
excluded by the Z° data flavour independent analyses below mh ~ 60 GeV. 
The region for mA < 10 GeV is difficult to exclude since the decay h°— >A°A° 
is usually dominant when kinematically allowed and the A*^ cannot decay to 
bb. Therefore the searches using b-tagging do not help in this region, while 
the sensitivity of the flavour independent h°Z° channels is too low to provide 
any exclusion. For a = and a — — 7r/4, the exclusion in this region is 
obtained by using the dedicated search for the process h°Z°— >A°A°Z° followed 
by A*'— s>cc, r"'"r~, gg and Z°^z/z/, e"^e~. By applying the Yukawa 
external constraint large tan/3 values are excluded. 

• The region 55 < mh < 60 GeV and mA > 75 GeV for a = in Figure 7(c) 
is not excluded for tan/3 < 1 since for the Z° data the h'^Z*^ cross-section be- 
comes too small to exclude it and most of the high energy flavour independent 
channels are only efficient for domains in which mh > 60 GeV. 

• The largest (mh, mA) excluded domain is for a = — 7r/4, where most of the 
parameter space points accessible at LEP are excluded, as can be seen in 
Figure 7(d). 

b) 0.4 < tan/3 < 1.0 (fighter grey area) and 1.0 < tan/3 < 40 (hatched area): 

• In Figure 7(c), as a consequence of the variation of the h°Z° production 
cross-section with tan/3, as discussed above, the mh < 10 GeV and 55 < 



mh < 60 GeV, ttt-a > 75 GeV regions are excluded for all values of in the 
tan/3 > 1.0 domain. For a = ih7r/2 the < 30 GeV region is excluded for 
all values of only in the tan/? < 1.0 domain. In the same figure, the region 
mh ~ 90 GeV and m\ > 60 GeV has an excess of data in the four-jet flavour 
independent h°Z° channel at 90 GeV and therefore cannot be excluded even 
for tan/3 > 1. 

• At a = and small values of tan/5 the production cross-section for the process 
e+e~—>h°Z° is highly suppressed. For mh > 60 GeV, constraining tan/3 > 1.0, 
larger excluded regions are obtained, as can be seen in Figure 7(c) (hatched 
area). In the same flgure, the unexcluded domain 88 < < 92 GeV, 
rrix > 60 GeV for tan/3 > 1 is due to the presence of candidates in the 
flavour independent four-jet channel in the year 2000 data [34]. 

• For a = ±7r/2. Figure 7(a), the region 100 < mh ^ 110 GeV and 35 < 

^ 50 GeV is unexcluded due to the presence of candidate events in the 
h^A'^— >-bbbb and h^A*^— s>bbr"'"r~ channels in the high energy data [35]. 

In Figure 8 the excluded regions in the (mh, mA_) plane independent of a are 
given together with the expected exclusion limits from MC studies. A particular 
(mh, niA) point is excluded at 95 % CL if it is excluded for all scanned values of 
tan/3 and a. Different domains of tan/3 are shown: 0.4 < tan/3 < 40 (darker 
grey area), 0.4 < tan/3 < 1.0 (lighter grey area) and 1.0 < tan/3 < 40 
(hatched area), for which enlarged excluded regions are obtained. A rectangular 
region 1 < mh ^ 55 GeV for 3 < ttia < 63 GeV is fully excluded at 95 % CL 
independent of a and tan/3. The scanned a > domain is new with respect to [6] 
and has the effect of restricting the exclusion for tan/3 < 1 to the kinematically 
accessible region for the h°A° production. 

For completeness, the excluded regions in the (mh, mA) plane for a < are 
given together with the calculated expected exclusion limits, is shown in Figure 9. 
The present study considerably extends the excluded (mh, mA) domain for negative 
values of a when compared with the study published in [6] . The previously excluded 
region of 1 < mh < 44 GeV and 12 < mA < 56 GeV is now enlarged to 1 < mh ^ 55 
GeV and 3 < mA ^ 63 GeV, for all tan/3 values for negative a. 

In Figures lO(a-d) the excluded regions in the (tan/3, mh) plane are shown for 
the chosen values of a, together with the calculated expected exclusion limits. A 
particular (mh, tan/3, a) point is excluded at 95 % CL if it is excluded for all scanned 
values of mA- There are two regions shown, the whole domain 3 GeV < mp^<2 TeV 
(darker grey area) and a restricted domain for which 3 < mA < 60 GeV (lighter 
grey area). The exclusion contours for mA < 60 GeV are larger for all a values, 
and entirely contain the 3 GeV < mA< 2 TeV excluded areas. 

Note that in Figure 10(b) the region tan/3 1 is unexcluded due to the suppres- 
sion of the h^'Z'' production cross-section as /3 — a ~ 0, while relatively low values 
of tan/3 can be excluded for 1 < mh < 60 GeV (darker grey area). Restricting the 
values of mA to be lower than 60 GeV improves the exclusion since the kinematical 
limit for h°A° production mechanism is never reached. 

In Figure 11 the excluded regions in the (mA, tan/3) plane are shown for different 
values of «, together with the expected exclusion limits. A particular (mA, tan/3, 
a) point is excluded at 95 % CL if it is excluded for all scanned values of mh- There 
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are three regions shown, corresponding to different mh domains that are subsets of 
one another, namely: 1 < mh < 110 GcV (darker grey area), 1 < mh < 90 GeV 
(hghter grey area) and 1 < mh < 75 GeV (hatched area). The lower the mh upper 
value analysed, the larger the excluded (mA, tan/3) region. 

In Figures 11 (a), (b) and (c) for niA < 10 GeV and tan/3 > 10 the excluded 
domains are a direct application of the Yukawa external constraint. 

Figures 11 (a) and (b) and Figures 10 (a) and (b) show a similar behaviour 
in the tan/3 excluded domains: for a = ±7r/2 {a = 7r/4) mostly low (large) tan/3 
values are excluded. 

In Figure 11 (c) for 1 < mh < 90 GeV the area with m\ > 45 GeV and tan/3 > 1 
is excluded as expected from the tan/3 > 1 contour in Figure 7 (c). In Figure 11 (d) 
the unexcludcd region in the dark grey contour for small tan/3 and tua < 50 GeV 
corresponds to the unexcluded region under the semidiagonal in Figure 7 (d). As 
can be seen in Figure 11 (d) for 1 < mh < 90 GeV, the area 24 < tua < 32 GeV 
for 3 < tan/3 < 14 is not excluded, but it is excluded when 1 < mh < 75 GeV as 
inferred from Figure 7 (d). 

In Figures 11 (b) and (d) for 1 < mh < 110 GeV the area 10 < uia < 25 GcV 
for tan/3 > 20 is excluded thanks to the optimisation of the e"'"e~^h°A''^bbbb 
analysis in the same kinematical region [35]. 

8 Conclusions 

New limits on the h°A° pair- production cross-section are obtained by the application 

of a new flavour independent h^A^ analysis at the highest LEP energies. A lower 
bound at 95 % CL is extracted along the diagonal at mh ~ mA ~ 71 GeV for 
— 1 assuming 100 % hadronic decays. The limit obtained by using the b-tagging 
analysis [35] and assuming 100 % decays into b-quarks, is at mh ~ mA ~ 81 GeV 
for = 1. 

A general analysis of the 2HDM(II) with no CP violation and no extra particles 
besides those of the SM and the five Higgs bosons has been performed using the 
Z*^, 183 and 189 GeV data together with the high energy data taken by OPAL in 
the years 1999 and 2000 at y/s = 192 — 209 GeV. Large areas of the parameter 
space of the model have been scanned. In the scanning procedure the dependence 
of the production cross-sections and branching ratios on the angles a and /3, calcu- 
lated with next-to-next-to-leading order QCD corrections and next-to-leading order 
electroweak corrections, has been considered. 

In addition to the standard OPAL b-tagging analyses, fiavour independent chan- 
nels for both the Higgs-strahlung process, e+e~— >h°Z°, and the pair-production 
process, e"'"e~— >h°A°, have been analysed, providing access to regions of parameter 
space in the 2HDM(II) where h° and A'' are expected to decay predominantly into 
up-type light quarks and gluons {e.g. a ^ 0). 

The 2HDM(II) parameter space scan, for 1 < mh < 130 GeV, 3 GeV < mA< 2 
TeV, — 7r/2 < a < 7r/2 and 0.4 < tan/3 < 40, leads to large regions being excluded 
at the 95 % CL in the (mh, mA) plane as well as in the (mh, tan/^) and (mA, tan/3) 
projections. The region 1 < mh < 55 GeV and 3 < mA < 63 GeV is excluded 
at 95 % CL both when restricting a < 0, as in a MSSM-like scenario, and in the 
general 2HDM(II) case, independently of a within the scanned parameter space. 
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The results obtained by adding the data taken by OPAL in the years 1999 and 
2000 substantially enlarge the excluded domains in the 2HDM(II) parameter space 
accessible by LEP. 
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Figure 1: The Monte Carlo mass grid on which the signal was generated for the 
year 2000 analysis. Similar grids are used for the four 1999 samples. 



OPAL 



2 5 

log(Ax (mj,,m^)) 



0.5 

lcos6, 



1 1000 

(U 

W 750 










- (c) 






• 


OPAL Data : 




□ 


Four-fermion ^ 




"I □ 


Two-fermion I 






10 pb Signal J 






m^=50 GeV : 






m^=60 GeV - 






(gggg) - 









0.2 



thrust 




-5 
log(Wcco3) 



-2 

l«g(WQCD) 



Figure 2: Distributions of the likelihood input variables for preselected events in the 
flavour-independent hadronic h°A° search: (a) log(Ax^(mh, wia)), (b) | cos 6'thrust I , 
(c) the event aplanarity A, (d) log(?/34), (e) J^, (f) (E^ax - ^mm)/v^, (g) the 
jet-charge-signed cos6'w, (h) log(Wcco3)) and (i) log(WQCD)- The dark-shaded 
histograms show the two-fermion (qq) background, the light-shaded histograms 
show the four-fermion background. The backgrounds are shown stacked atop one 
another. The dashed histograms show a 10 pb signal of e+e~— s>h°A''^ gggg with 
mh = 50 GeV and niA = 60 GeV, and the points with error bars show the OPAL 
data collected in 1999 and 2000. 



OPAL 



^800 



600 
400 
200 




1 800 
W600 



400 



(a) 



m^=50GeV m;^=100GeV 



10 pb Signal- 



m^=50 GeV m^=100 GeV ^^^y 



• OPAL Data 
^ □ Two-fermion 
I I Four-fermion 



(bbbb) : W 




10 pb Signal- 
(gggg) 




log(Ax (mj^,m^)) 



0.4 0.6 
Likelihood 



log(Ax (m. ,m,)) 



0.2 0.4 0.6 
Likelihood 



log(Ax (m. ,m )) 



0.4 0.6 
Likelihood 



m^=5b GeV m^=60 GeV 





Figure 3: Distributions of the log(Ax^(mh, mA)) variable and the likelihood out- 
put variable for three test-mass hypotheses, (mh, mA) = (50 GeV, 100 GeV), 
(50 GeV, 60 GeV), and (30 GeV, 60 GeV). OPAL data are shown with the points, 
light-shaded histograms show the four-fermion Standard Model background expec- 
tations, and dark-shaded histograms show the two-fermion (qq) Standard Model 
background expectations. The two background rates are shown stacked in all his- 
tograms. In the likelihood histograms, a 1 pb signal assumed to decay to gggg is 
shown with the hatched histograms, added on top of the background sum. For the 
log(Ax^(?7ih, "h^a)) distributions, a 10 pb signal assumed to decay to bbbb is shown 
with dot-dashed histograms which are not added to the background, while a 10 pb 
signal assumed to decay to gggg is shown with dotted histograms, also not added 
to the background. 
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Figure 4: The value of 1 — CLb as a function of (mt, mp^). No model has 1 — CLb < 
0.01, and one or more excesses with significance 1 — CLb < 0.01 are expected due 
to the large number of independent searches conducted. 
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Figure 5: The median expected limit on the couphng c^. Both the h° and the A° 
are assumed to decay hadronically 100% of the time. 
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Figure 6: The limit on the couphng c^. Both the h° and the A° are assumed to 
decay hadronically 100% of the time. 
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Figure 7: Excluded regions in the (mh, mx) plane, (a)-(d), for a = ±vr/2, 7r/4, 
and — vr/4, respectively, together with the expected exclusion limits. A particular 
(mh, rriA, a) point is excluded at 95% CL if it is excluded for all scanned values of 
tan/3. Three different domains of tan/3 are shown: the darker grey region is excluded 
for all values 0.4 < tan/?< 40; additional enlarged excluded regions are obtained 
constraining 0.4 < tan/? < 1.0 (lighter grey area) or 1.0 < tan/?< 40 (hatched area). 
Expected exclusion limits for 0.4 < tan/?< 40 are shown as a dashed line. 
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Figure 8: Excluded (mA, region independent of a, together with the expected 
exclusion limit. A particular (mA, point is excluded at 95% CL if it is excluded 
for 0.4 < tan/? < 40 (darker grey region), 0.4 < tan/5 < 1.0 (lighter grey region) 
and 1.0 < tan/3 < 40 (hatched region) for any a. Expected exclusion limits for 0.4 
< tan/3< 40 are shown dashed line. 
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Figure 9: Excluded (mA, rrih) region for — 7r/2 < a < 0, together with the 
expected exclusion limit. A particular {rrij^, mh) point is excluded at 95% CL if it 
is excluded for 0.4 < tan/5 < 40 (darker grey region), 0.4 < tan/? < 1.0 (lighter grey 
region) and 1.0 < tan/3 < 40 (hatched region). Expected exclusion limits for 0.4 < 
tan/?< 40 are shown dashed line. 
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Figure 10: Excluded regions in the (tan/?, mh) plane, (a)-(d), for a= ±vr/2, 7r/4, 
and — vr/4, respectively, together with the expected exclusion limits. A particular 
(mil, tan/?, a) point is excluded at 95% CL if it is excluded for all scanned values of 
m-A- The two regions shown correspond to the whole domain 3 GeV < mA< 2 TeV 
(darker grey area) and a restricted domain for which 3 < rriA < 60 GeV (lighter 
grey area). The exclusion regions for mA < 60 GeV entirely contain the 3 GeV 
< mA< 2 TeV excluded areas. Expected exclusion limits are shown for 3 GeV < 
"^A< 2 TeV (dashed line). 
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Figure 11: Excluded regions in the (mA, tan/?) plane, (a)-(d), for a= ±7r/2, 7r/4, 
and — 7r/4, respectively, together with the calculated expected exclusion limits. A 
particular (mA, tan/3, a) point is excluded at 95% CL if it is excluded for all scanned 
values of mh- The three contours correspond to 1 < mh < 110 GeV (darker grey 
area) and 1 < mh < 90 GeV (lighter grey area). Expected exclusion limits are 
shown for 1 < mh < 110 GeV (dashed line). 
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